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ABSTRACT 

We  discuss  the  use  of  homogenization  techniques  to  derive  approximate  models  with  sim¬ 
ple  geometry  for  physical  models  of  grids  and  trusses  which  have  a  complex  geometry  that 
gives  rise  to  computational  difficulties.  Our  presentation  is  in  the  context  of  inverse  or  pa¬ 
rameter  estimation  problems  for  composite  material  structures  with  unknown  characteristics 
such  as  stiffness  and  internal  damping.  We  present  the  necessary  theoretical  foundations  for 
this  approach  and  discuss  comparison  of  modal  properties  of  the  resulting  homogenization 
model  for  a  two-dimensional  grid  structure  with  modal  properties  observed  in  experiments 
with  this  grid. 
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1.  Introduction 


We  report  here  on  a  part  of  our  continuing  efforts  on  the  development 
of  high  fidelity  dynamic  models  for  composite  material  structures.  The  fo¬ 
cus  of  our  investigations  has  been  on  models  to  be  used  in  estimation  and 
control  of  large  flexible  structures  mainly  intended  for  use  in  space  (e.g.,  an¬ 
tennas,  platforms,  solar  panels,  experimental  arrays,  etc.).  At  present  there 
is  a  reasonably  adequate  understanding  of  the  dynamics  of  beams  and  plates 
made  from  known  materials  such  as  aluminum  alloys.  Our  recent  efforts  (see 
[Bln]  and  the  references  therein)  involved  models  and  methods  to  determine 
material  parameters  in  composite  material  structures  with  simple  geometry 
(beams  with  attached  solid  bodies  or  solid  plates).  These  result  in  inverse 
problems  that  are  computationally  tractable  as  long  as  the  physical  geometry 
is  relatively  simple.  However,  substantial  difficulties  arise  in  cases  involving 
more  complex  geometries  such  as  grids  (which  may  be  viewed  as  plates  with 
many  holes)  and  trusses  (solid  columns  from  which  most  of  the  material  is 
removed  in  some  periodic,  regular  fashion).  In  these  cases  the  difficulties 
associated  with  unknown  composite  material  characteristics  such  as  stiffness 
and  internal  damping  are  combined  with  severe  difficulties  related  to  com¬ 
putational  grid  selection  for  a  domain  that  is  mostly  holes  or  perforations. 

The  purpose  of  this  paper  is  to  outline  one  possible  methodology  for 
dealing  with  these  structures  of  complex  geometry  and  to  report  on  some 
of  our  initial  investigations  in  this  regard.  This  methodology  is  based  on 
ideas  from  homogenization  and  requires  that  the  structures  be  highly  periodic 
(many  perforations  repeated  in  a  regular  pattern)  and  sparse  in  material  (i.e., 
a  grid  or  truss  with  thin  members). 

We  formulate  our  ideas  in  the  context  of  an  inverse  problem  approxi¬ 
mation  framework  that  has  been  developed  for  problems  involving  simple 
geometries.  A  summary  of  this  theoretical  framework  is  presented  in  Section 
2.  In  Section  3  we  describe  a  particular  grid  structure  which  we  have  used  in 
our  experiments  and  indicate  how  to  formulate  a  direct  physical  model  for 
this  grid  within  the  framework  of  Section  2.  This  model  leads  to  computa¬ 
tional  methods  that  are  inherently  intractable  for  estimation  of  parameters 
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(and  also  for  the  ultimately  desired  control  investigations)  due  to  the  under¬ 
lying  computational  domain. 

In  Section  4  we  present  a  summary  of  results  for  a  homogenization  proce¬ 
dure  that  approximates  (in  a  nonstandard  way)  the  original  direct  physical 
model  on  a  perforated  domain  by  a  homogenized  model  on  a  domain  that  is 
very  simple  (the  perforated  domain  with  all  of  the  perforations  filled  in).  The 
resulting  approximate  model  is  also  in  a  form  to  fit  into  the  inverse  problem 
framework  of  Section  2.  Finally,  in  Section  5  we  report  on  our  initial  efforts 
on  validation  of  this  approximate  model  by  comparing  experimentally  ob¬ 
served  modal  properties  of  the  grid  with  those  possessed  by  the  homogenized 
model. 


2.  Review  of  Theory  for  Second  Order  Systems 


In  this  section  we  give  a  brief  summary  of  the  theoretical  background 
necessary  for  a  rigorous  discussion  of  estimation  problems  for  second  order 
systems.  Detailed  discussions  can  be  found  in  [BI],  [BR 1] ,  [BR2],  [BW],  [Bj, 
[BK], 

Let  V  and  H  be  complex  Hilbert  spaces  with  V  continuously  and  densely 
embedded  in  H.  We  may  then  formulate  a  Gelfand  triple  V  c— ♦  //  >  V*  with 

duality  pairing  (  ,  )v\v  (e.g.,  see  [W]),  which  we  shall  denote  by  (  .  )  when 
no  confusion  results.  Let  ( Q,d )  be  a  compact  metric  space  of  admissible 
parameters  </.  We  consider  the  parameter  dependent  abstract  second  order 
system 


(2.1)  u(t)  +  A2(q)u(t)  -f  A\(q)u(t)  —  f{t,q)  t>  0 

u(0)  =  u0,  «(0)  =  t’o, 

where  the  operators  A{(q)  G  C{V,Vm),i  =  1,2,  arise  from  parameter  de¬ 
pendent  sesquilinear  forms  cr,(q)  :  V  x  V  — ►  C  which  represent  generalized 
stiffness  (crx )  and  damping  (<72),  respectively.  More  precisely,  we  consider  the 
equation 

(2.2)  {u{t),<t>)  +  a2{q)(u{t),<p)  +  ai{q){u{t),<p)  = 

for  all  <t>  G  V,  where  C\  is  symmetric  and  oq.oy  satisfy  the  following  conditions 
(a  theory  for  damping  forms  ay  that  satisfy  weaker  conditions  than  (B). 
namely  //-semiellipticity,  can  be  found  in  [HI]): 

(A)  Parameter  Continuity:  For  each  0  G  V  there  is  a  function 

c  (*,  •;  <f)  on  Q  x  Q  satisfying  ;(q,  p\  q)  — *  0  as  q  — >  p  in  Q  such  thci 
for  each  if  G  V  we  have 

-  <7(p)(o,  t’)|  <  s(y,p;<?)K’|v. 

(B)  K-ellipticity:  There  exists  ci  >  0  such  that  for  all  q  G  Q  and  all 
<t>  G  V  we  have 

Be  er(e/)(0.  o)  >  C||o|f.. 


(C)  Boundedness:  There  exists  c2  >  0  such  that  for  all  q  £  Q  and  all 
<j>,  if  6  V  we  have 


Under  these  conditions  on  oq  and  <x2,  there  exist  operators  A,(q)  6 
C(V,  V*)  such  that 

°i{q){<t>,  V’)  =  { Ai(q)<l> ,  0)  ^,0  6  V. 

Hence  equation  (2.1),  or  equivalently,  (2.2),  is  to  be  interpreted  as  an  equation 
in  V*.  Following  standard  practice,  we  may  rewrite  (2.1)  in  first  order  vector 
form  in  the  coordinates  w  =  ( u,ii )  and  use  semigroup  considerations  when 
discussing  solutions.  To  this  end,  we  define  Tt  =  V  x  H  and  V  -  U  x  V 
and  note  that  V  <->  Tt  »  V'  also  forms  a  Gelfand  triple.  We  also  define  the 
sesquilinear  form  <r(q)  :  V  x  V  — >  C  by 

(2.3)  <r(q)((t,tf)>(q,0)  =  ~(4’,h)v  +  <7, (?)(<£>£)  +  ^(tfXVbO 

so  that  (2.2)  may  be  written 

(*»(<)»  x>  +  <r(?)(M0>x)  =  {F{t,q),x) 


(2.4) 


w(0)  =  (u0,  v0) 


for  w(t)  =  (u(f),u(f))  and  x  =  {4>,  i>)  in  V  with  F(f,q)  =  (0,  f{t,q)).  Equiv¬ 
alently,  we  may  write  the  equation  in  V*  as 


(2.5)  w(t)  =  A{q)w(t)  +  F{t,q) 

where  cr(q)(x,  C)  =  (-^(tfhsO  with 


-4(<?)  = 


0 


/ 


l  -am  -am 

Since  <7i  and  a2  satisfy  (B)  arid  (('),  it  is  readily  established  that  A(q)  is 
the  infinitesimal  generator  of  an  analytic  semigroup  T(t\q)  on  V.7d  and  V". 


1 


Moreover,  the  unique  solution  w  £  L2(0,T;hi)  of  (2.5)  for  w0  =  (n0,u0)  £ 
hi  =  V  x  H  and  /  £  Z/2(0,T;  7/)  is  given  by 

(2.6)  w(t\ q)  =  T(t;  q)w0  +  f  T (t  -  s\q)F(s,q)ds. 

Jo 

For  computational  purposes,  one  must  consider  approximation  schemes 
for  (2.6).  Let  H*  be  a  family  of  finite  dimensional  subspaces  of  H  satisfying 
H*  C  V  and  the  condition 

(Cl)  For  each  (f>  6  V,  there  exists  <}>*  £  H*  such  that  \<j>  —  q>*  (v  — +  0  as  N  — *  oo. 

Let  hi*  =  H*  x  HN  and  let  PN  be  the  orthogonal  projection  of  hi  onto  hi* . 

Moreover,  let  AN(q)  €  C{hi *)  denote  the  operator  obtained  by  restricting 
a(q)  of  (2.3)  to  hi*  x  hi* \  i.e.,  for  xsC  €  hi* , 

^(?)(X.C)  =  (~A*(q)x,  ()■ 

We  denote  by  T*(t;q)  =  eA  ^  the  corresponding  analytic  semigroups  on 
hi*.  The  approximating  systems  for  (2.6)  are  then  defined  by 

(2.7)  w*(t ;  q)  =  TN(t;  q)P* Wq  +  /'  T*(t  -  s;  q)PN F(s,  q)ds. 

Jo 

From  an  application  of  the  Trotter- Kato  theorem  on  convergence  of  semi¬ 
groups,  one  can  readily  obtain  that  under  conditions  (Cl),  (A),  (B),  (C). 
we  have  for  each  \  G  hi,T*  (t;  q*)P* \  — >  T(t-,q)x  in  hi,  uniformly  in  t  on 
compact  intervals,  whenever  {q* }  is  an  arbitrarily  chosen  sequence  in  Q  with 
q*  — >  q.  Moreover,  using  the  theory  of  analytic  semigroups,  one  can  argue 
that  for  each  \  £  hi  and  each  positive  integer  k,  A* (q* )kT*  (t;  q* ) P*  \  — > 

A(q)kT(t;  q)x  in  hi,  again  uniformly  in  t  on  compact  intervals. 

It  follows  immediately  that  w*(t;q*)  of  (2.7)  converges  in  hi  to  ir(t;q)  of 
(2.6),  with  the  convergence  being  uniform  in  t  on  compact  intervals.  Further¬ 
more,  for  /  sufficiently  smooth  (e.g.,  /  £  CM((0,  71],  //)  suffices),  we  find  that 
w*(t]q *)  — >  w(t;q)  in  hi  as  N  —>  oc,  uniformly  in  t  on  compact  intervals. 

Thus,  we  see  that  the  first  component  u*{t]q*)  converges  to  ut(t:q )  in  1 
while  the  approximate  acceleration  u*t(t:q*)  converges  in  II  to  utt[t\q).  With 
more  smoothness  on  /  (e.g.,  /  £  C2( [0.  7’].  II)).  we  obtain  that  u*  (/;  </x )  con¬ 
verges  in  V  to  att(t\q)  as  ,Y  — ►  oo. 


If  V  embeds  continuously  in  C(Q)  (this  is  the  case,  for  example,  if 

V  C  H2(Vl)  where  f l  C  R 1  or  0  C  R 2),  then  the  above  results  lead  to 
pointwise  convergence  (in  the  spatial  variables  as  well  as  time)  of  approxi¬ 
mations  to  the  displacement,  velocity  or  acceleration  in  dynamic  problems 
involving  structures  such  as  beams  and  plates.  These  approximation  results 
are  precisely  those  needed  to  treat  certain  questions  arising  in  inverse  or 
parameter  estimation  problems  for  such  structures. 

To  be  more  specific,  suppose  we  are  given  a  set  of  observations  in 
the  observation  space  Z  along  with  an  observation  map  F(</)  from  T>  C 
C([0,T],H)  to  Z  for  the  system  (2.1).  Thus,  2  is  a  set  of  observations  or 
data  for  Tu{q)  where  u(q)  is  the  solution  to  (2.1).  A  least  squares  estimation 
problem  then  consists  of  finding  q  €  Q  which  minimizes  over  Q  the  criterion 

(2.8)  4>(g)  =  M9)-2||, 

subject  to  (2.1).  Typical  examples  included  in  this  setting  are  problems 
for  beams  and  plates  where  acceleration  measurements  (accelerometers)  are 
available.  In  the  case  of  a  beam  of  length  a,  we  might  have  observations  zjk 
of  the  acceleration  utt{tj,xk-,q)  at  several  times  tj,j  =  1,2, in  [0.  T] 
and  at  several  locations  xk  €  (0, =  1,2 along  the  beam.  Then 

V  C  /f2(0,a)  in  the  usual  formulation  and  V  =  C2([0,71],  V'),  Z  =  RJh  with 

r  u(q)  =  { utt(tj,Xk‘,q )  :  j  =  k  =  1,  —  ,  AT}  so  that 

$(?)  =  H\ uu{tj,xk\q)  -  zjk\2. 

For  a  rectangular  plate  occupying  the  region  (x,  y)  €  12  =  (0,  a)  x  (0,  b).  use  of 
accelerometers  at  locations  (x* .,  yk),  A;  =  1, . . . ,  I\,  would  lead  to  observations 
{zjk}  for  {uu(tj,xk,yk;q)}.  With  V  C  //2( Q),V  =  C2(\ 0,7'], V)  and  2  = 
RJh  we  have  T u(q)  =  {uu{tj,xk,  yk\q)  :  j  =  1,  - .  - , k  =  1, .  •  • ,  A , }  and 

$(<7)  =  H  Wtt((j,xk,yk:q)  -  zik |2. 

j-k 

The  corresponding  approximate  problems  are  formulated  using  the  crite¬ 
rion 

(2.9)  <t>X(q)  =  |IVv(c/)  - 


() 


for  solutions  w^(q)  =  (uN (q) ,  iiN {q))  of  (2.7).  One  then  seeks  a  parameter 
qN  €  Q  that  minimizes  $N(q)  over  q  6  Q  subject  to  (2.7). 

Among  the  important  questions  related  to  such  approximate  problems 
are  those  pertaining  to  parameter  convergence ;  i.e.,  for  a  given  fixed  set  of 
data  z  £  2,  do  optimal  parameters  qN  for  (2.9)  converge  in  some  sense  to 
an  optimal  parameter  for  (2.8)?  More  generally,  one  can  also  incorporate 
continuous  dependence  of  the  optimal  estimates  on  the  data  by  employing 
the  concept  of  method  stability.  If  one  is  given  a  sequence  {zm}  of  data 
sets  that  converge  in  Z  to  a  data  set  2°,  and  one  denotes  by  qh(zm)  and 
q(z° )  the  optimal  parameters  (in  general,  these  are  sets)  for  (2.9)  and  (2.8) 
corresponding  to  zm  and  z°,  respectively,  then  method  stability  requires  that 
qN(zm)  converges  (in  some  appropriate  setwise  sense)  to  q(z°)  as  N,  m  — >  oo. 
These  issues  are  carefully  discussed  in  [B],  [BK],  where  it  is  shown  that  to 
insure  both  parameter  convergence  and  method  stability,  it  suffices  to  argue 
that  T uN(qN)  — >  ru(<7)  in  Z  for  any  sequence  {<7^}  in  Q  with  qN  — >  q.  From 
the  discussions  above  one  thus  has  that  the  conditions  (A),  (B),  (C)  on 
and  <7 2  and  (Cl)  along  with  a  compact  parameter  space  ( Q,d )  are  sufficient 
to  treat  these  questions  in  the  case  of  accelerometer  data  for  beams  and 
plates. 
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3.  The  AFAL  Grid  Structure 


We  consider  a  rectangular  plate  perforated  with  rectangular  holes  as  de¬ 
picted  in  Figure  3.1.  A  similar  plate  has  been  the  subject  of  numerous 
experimental  investigations  at  the  Air  Force  Astronautics  Lab  (now  a  part  of 
the  Phillips  Lab)  at  Edwards  Air  Force  Base.  Using  the  Love-Kirchhoff  plate 
theory  (the  2-dimensional  analogue  of  the  Euler- Bernoulli  theory  for  beams), 
one  can  model  such  a  grid  structure  as  a  second  order  system  of  the  form 
(2.1)  so  that  conditions  (A),(S),(C)  of  the  previous  section  are  satisfied. 
We  summarize  here  previous  findings  and  refer  the  reader  to  [BR1],  [BR2], 
[Rj  and  the  references  therein  for  more  detailed  discussions. 


Figure  3.1 

We  assume  that  the  plate  is  cantilevered  along  the  a-i  axis  (at  .r2  =  0)  as 
depicted  in  Figure  3.1  and  is  free  on  the  other  edges  as  well  as  the  edges  of  the 
holes  (perforations).  If  we  use  the  Love-Kirchhoff  theory  and  assi.no  Kelvin- 
Voigt  damping  as  well  as  viscous  (air)  damping,  the  transverse  displacements 
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u(t,x  1,^2)  at  time  t  and  location  x  =  (x^a^)  in  the  perforated  domain 
f2per  =  (0,  a)  X  (0,  b)  —  {  perforations  }  satisfy  a  second  order  (in  time)  system: 


d2u  du  d2Ml  nd2M12  d2M2 

(3.i)  ph—+i—  +  +  + 


—  /  t  >  0,  (xj ,  x-2)  €  0 


peri 


where  p  is  the  mass  density,  h  is  the  thickness  of  the  plate,  and  A/1,  A/2,  and 
A/12  are  the  bending  moment  about  the  x'i  axis,  the  bending  moment  about 
the  X2  axis  and  the  twisting  momcn*  respectively.  If  we  include  Kelvin- Voigt 
damping  in  these  moments,  we  find 

El  ( d2u  #2u1  C£>I  (  d3u  d3u  ) 

1  —  i/2  |  dx\  dx\  J  +  1  —  v2  {  dx\dt  ^  dx2dt  J 

El  ( d2u  d2u  \  cDI  f  d3u  d3a 

1  —  u2  {  dx\  ^  dx\  J  ^  1  —  i/2  \  dx2dt  dx\dt 

El  02u  cd  I  d3u 
- - -  — - 

1  +  v  dx\dx2  1  +  v  dxidx2dt 

where  El  and  col  are  the  usual  stiffness  and  Kelvin- Voigt  damping  coeffi¬ 
cients  and  v  is  Poisson’s  ratio.  The  plate  is  clamped  along  the  Xi  axis  where 
we  must  have  the  (essential)  boundary  conditions 

du 

(3.3)  u  —  — —  =  0,  for  x-2  =  0,  0  <  Xj  <  a,  t  >  0. 

dx  2 

Since  the  plate  is  free  on  its  other  edges  and  the  hole  edges,  it  must  satisfy  zero 
moment  and  zero  shear  conditions  on  these  edges.  These  result  in  natural 
boundary  conditions  given  by 


A/1  = 

(3.2)  M2  = 

A/12  = 


(3.4) 


Ml  =  0, 


DM1  93A/12 
(hr  1  ~  dx-2 


on  edges  parallel  to  the  X2  axis,  and 


(3.5) 


A/2  =  0, 


OM2  ,/h\/12 

Ox  2  ~  Ox  1 


on  edges  parallel  to  the  xj  axis. 


The  system  (3.1)  -  (3.5)  can  be  readily  formulated  in  the  sesquilinear 
framework  of  Section  2  so  that  the  results  given  there  are  applicable.  To 
that  end,  define 


Vper  =  Ht{nper)  =  {(pe  H2(Q per)  :  <j>  =  - -  =  0  On  X2  =  0} 

OX  2 


and 


//per  —  L2  (flper). 


Then  lper  //per  c— ♦  V p’er  forms  a  Geifand  triple.  Moreover,  we  may  define 
sesquilinear  forms  crfer(q)  :  lper  x  Vper  — +  C,  i  =  1,2  by 


per 


(9)(0>0)  = 


(3.6) 


{9i(0,ll  0,11  T'f',2'2  0,22  )  +  9i93(0,l  1  0,22  +  0,22  0,11  ) 


+29,(1  -  <73)0,12  0,i2  }0G 

where  9  =  (91,92,93,94)  =  (7— 7.T^4,b7)  and  0,tJ  =  It  is  not 

1  —  1  12  *  OXiOXj 

difficult  to  verify  that  crf<r(q)  satisfies  conditions  (/l),  (/?),  (C)  whenever  Q  is 
a  compact  subset  of  C(flper, /t"1)  H  {9  :  9,  >  c  >  0,i  =  1,2, 3;  94  >  0, 93  <  1}. 
The  system  may  then  be  written  for  u(t)  G  lper  as 

(3.7)  (p/a/it,i/’)  +  cr(jer(9)(u(O,  </’)  +  crrr(9)(«(O,0)  +  (94 «(Oi0)  =  </,  0) 

for  all  0  G  V0,r.  Given  initial  data  u(0)  =  wo,u(0)  =  i'o  with  n0  G  Vpcr,  vq  G 
// ptr,  the  methods  outlined  in  Sectie  2  can  lie  used  to  develop  direct  numer¬ 
ical  methods  as  discussed  in  detail  in  [BR2].  While  the  theoretical  aspects 
of  approximation  in  parameter  estimation  problems  are  relatively  straight¬ 
forward  conceptually  given  the  theory  of  Section  2,  the  problems  are  com¬ 
putationally  quite  demanding.  This  is  true  even  though  the  computational 
domain  flper  is  among  the  simplest  and  most  regular  of  those  arising  in  grid 
and  truss  structures. 

The  computational  difficult  k\:  alluded  to  here  motivate  one  to  develop 
alternative  methods  for  modeling  of  grid  and  truss  structures  for  the  purposes 
of  identification  and  control. 
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4.  Homogenization  for  Grid  Structures 


Given  the  difficulties  associated  with  the  model  for  the  grid  described  in 
the  previous  section,  it  is  of  interest  to  approximate  grid-like  and  truss-like 
structures  with  a  model  that  retains  physical  fidelity  but  is  computationally 
more  tractable  than  a  direct  physical  model.  One  such  approximate  model 
can  be  formulated  using  homogenization  techniques  which  were  original! v 
developed  to  model  composite  material  structures.  Here  we  outline  results 
for  a  grid  structure  similar  to  that  of  Section  3.  Details  for  this  grid  and 
other  lattice/truss  derivations  can  be  found  in  [BCR],  [R],  [BLR],  [CS].  [CD], 
and  [SR]. 

The  basic  idea  in  our  use  of  homogenization  is  to  approximate  a  model 
such  as  (3.1)  -  (3.5)  or,  equivalently,  (3.7).  by  a  homogenized  model  (HI'.  ';  on 
a  domain  in  which  the  "holes"  or  perforations  have  been  "filled  in".  That 
is,  instead  of  a  perforated  domain  D,)cr  as  described  in  Section  3.  we  wish  to 
compute  on  the  domain  Q  =  (0,a)  x  (0.6).  The  model  associated  with  this 
domain,  even  if  it  is  plate-like  in  form.  will,  of  course,  be  nonphysical.  The 
coefficients  in  the  distributed  system  will  be  nonphysical  parameters. 

The  homogenization  procedures  we  outline  here  require  periodicity  in  the 
structure  (many  regularly  placed  "holes")  as  well  as  sparseness  of  material 
in  the  structure.  Thus  we  consider  a  grid  structure  like  that  in  figure  3.1. 
but  with  manv  "holes"  or  bavs.  We 


Xi 


x2 

— iif  ”  "T. 


I  I 

I  I 


X  =  (xi,x2) 


y  2 


Y  =  [0,£i]  X  [0,e2] 
=  K*UT 


Figure  4.1 


as 

consider  a  typical  cell  as  depicted  in  Figure  4.1.  We  map  this  typical  cell  to 
a  fixed  size  cell  Y  =  [0,<d]  x  [0,  £2]  via  a  mapping  y  —  x/e  so  that  the  original 
cell  has  dimension  et\  by  et2.  As  £  -+  0,  periodicity  in  the  original  fixed 
domain  flper  increases.  As  /z  — *  0,  the  thickness  of  the  members  in  the  grid 
decreases.  Our  approximate  equation  on  Q  will  be  obtained  by  taking  limits 
as  e  -+  0,/<  — *  0.  Of  course,  the  approximate  limit  model  can  be  expected  to 
be  a  better  approximation  to  the  actual  grid  structure  if  the  grid  has  a  large 
number  of  bays  and  thin  members. 

To  facilitate  our  discussions,  we  parameterize  the  grid  structure,  denoting 
the  perforated  region  by  0,£^,  and  consider  a  model  of  which  (3.7)  is  a  special 
case.  (For  convenience,  we  drop  the  viscous  damping  term  for  the  discussions 
in  this  section.)  We  define  H£4l  =  L2{Qc,n)  afid  V£tP  =  //fc2(fif,(1)  =  {o  € 
=  0  on  x2  —  0}.  The  generalized  stiffness  and  damping 
sesquilinear  forms  are  defined  as 
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for  <f>,  ip  6  Ve,n-  Here  and  throughout  we  adopt  the  notation  4>uj  f°r  ofdx  ■ 
We  assume  that  the  functions  defined  by 

aljkh  (jJ  =  aijidx)  and  bi]kh  ^  =  b]?kh{x)  for  x  <E 
are  K-periodic  and  that  initial  data 

are  given.  We  then  consider  the  system  for  ue'u(t)  €  satisfying  for 

0  €  Vtllt 

(4.3) 

{phu\f(t),  xp)  +  t  6  (0,T), 

u'^O)  =  u^,  uM  0)  =  vS* 

We  wish  to  take  limits  as  e  — >  Q,p  — ►  0  in  this  system.  Note  how¬ 
ever  that  as  e  and  p  change  so  does  the  domain  hence  limit  taking 

is  not  a  straightforward  process.  We  need  to  first  extend  functions  such  as 
ue,tJ-(t),  /,  etc.  to  the  entire  domain  Q.  We  can  do  this  by  assigning  the 
value  zero  to  functions  in  the  region  Q  —  We  shall  denote  by  g  this  ex¬ 
tension  by  zero  of  any  function  g  £  This  results  in  g  6  L2(H).  Such 

an  extension  is  not  so  useful  if  we  wish  to  deal  with  functions  g  6  //2(fL<M) 
(such  as  u£,^(t))  for  which  the  H2  smoothness  is  to  be  preserved.  For  such 
an  extension,  special  extension  operators  must  be  formulated. 

We  shall  proceed  in  two  steps.  We  fix  p  at  some  value  and  take  limits 
as  e  — ♦  0,  obtaining  an  intermediate  “homogenized”  solution  u M  satisfying  a 
certain  intermediate  equation.  We  then  take  the  limit  as  //  — >  0,  obtaining 
the  desired  homogenized  system  (HE). 

Let  p  >  0  be  fixed  and  let  ue  =  ue’11  be  the  solution  to  (4.3)  where  we 
shall  in  our  notation  temporarily  suppress  the  dependence  on  the  fixed  value 
of  p.  Then  we  have 

Lemma  4.1.  There  exists  an  extension  operator  Pc 

P'  e  C(L°°(0,  T ;  Vf,„),  Loc(0,  T ;  H^(ft))) 


13 


and  a  function  u  =  u M  suck  that  for  some  sequence  en 


0  we  have 


P'nuen  — >  u  in  I°°(0,T;  //b2(D))  weafc  * 

( Penuen)t  =  Penuct”  ->  it*  in  L°°(0,T;  L2(D))  weak  '. 


With  additional  assumptions  on  the  initial  data  in  (4.3),  we  can  determine 
a  limit  equation  that  it  =  uu  of  Lemma  4.1  satisfies.  Define  0  =  |5  |/|1  |, 
where  |F|  denotes  the  measure  (area)  of  Y.  Then  we  have 

Theorem  4.1.  Let  uc  =  ii6’*  he  the  solution  of  (4.3)  and  suppose  that 

f  €  L2{ile^)  so  that  f  €  L2(Ll).  Moreover,  suppose  there  exist  u0  €  H^(Q) 

and  vo  €  L2(Q )  such  that  the  initial  data  for  (4-3)  satisfy 

Uo’M  — >  u0  as  e  — >  0,  weakly  in  L2(Ll) 

Vq"  ->  u0  as  e  —>  0,  weakly  in  L2(Q). 

Then  the  convergence  of  Lemma  4-1  holds  with  the  limit  function  u(i)  = 
u^(t)  e  Hl(Tl)  satisfying  the  homogenized  system 

[HE") 

(i 9phuu(t),if )  +  o(t)(u(-),rf)  =  (Of, if)  for  t  €  (0 ,T)  and  all  4’  €  //fc2(D) 

u(0)  =  uo/0,  iq(0)  =  vo/0, 

corresponding  to  the  homogenized  hysteresis  sesquilinear  form 

(4.4)  a(t)(u{-),4')  =  /  f  Pijkh{t  -  T)u,kh{T,x)if>,ij{x)dT<L r. 

The  coefficients  pijkh{t )  =  L~'  [ptjkh{s)}  arc  the  inverse  Laplace  transforms  of 
the  functions 

(4.5) 

Ptjkh(s)  =  L[pljkk}{s)  = 

■—  /  { aljkh{y )  +  sb,jkh{y)  -  YL(a(mkh{y)  +  sb(mkh{y))x,cm  (S-U )Wl/ 

D  I  4v*  <,m 
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where  the  functions  y'J(s,  y)  are  Y -periodic  solutions  of 

I  X  ( a(mkh{y )  +  sb{mkk(y))(x‘J(s,y)  -  V,J{y)),em  4’,kh  (y)dy  =  0 

JY  (,m,k,h 

for  all  Y -periodic  4’  €  H2(Y*).  Here  VlJ(y)  =  \ytyj. 

We  note  that  the  general  homogenized  plate  equation  (for  //  fixed)  for 
a  grid  with  damping  involves  a  time  hysteresis  functional  of  the  solution. 
It  is  instructive  to  consider  the  special  case  of  the  grid  with  no  damping; 
i.e.,  bi-jkh  =  0.  In  this  case  we  find  that  pijkh(s)  =  L^puR-s)  is  actually  a 
constant,  say  ptJkh(s)  =  alJkh  with 

(4.6)  aijkh  =  jpr  j  { atJkh(y )  -  X  a(mkh{y)x'em  ( y))dy 

*  '  £,m 

where  x u  is  the  solution  of 

[  X  aemkh(y){x,J(y)  -  E'J(y)),tm  4\kk  ( y)dy  =  o 

JY'  l,m,k,h 

for  all  T-periodic  xf  £  IP(Y“).  Recalling  that  the  inverse  Laplace  transform 
of  a  constant  is  that  constant  times  the  Dirac  delta  function  (i.e.,  L~l[at]kh\  = 
atjkhh),  we  find  the  sesquilinear  form  (4.4)  reduces  to 

a(t)(u(-),xp)  =  I  ]T  atJkhu,kh  (t,x)4>UJ  {x)dx, 

or 

(4.7)  X  dijkh<t>,kh  (•r)v/’,iJ  [x)dx. 

i,j,k,h 

Thus,  the  hysteresis  in  the  model  is  a  result  of  nontrivial  damping  terms  in 
the  grid  model.  Moreover,  we  point  out  that  the  coefficients  in  (4.7).  and 
hence  the  coefficients  in  (HE  tl )  for  grids  without  damping,  are  not  simple 
averages  of  the  original  coefficients  atJkh  of  (4.3).  Rather,  the  homogenized 
coefficients  are  averages  of  the  original  coefficients  over  the  structure  plus 
some  correction  terms  as  shown  in  (4.6). 

We  proceed  to  the  next  step  by  letting  p  — >  0  in  ( HE " )  and  the  associated 
equations  given  in  Theorem  4.1.  To  facilitate  our  discussions,  we  restrict 


our  considerations  to  the  constant  coefficient  case;  i.e.,  a>3kh  and  bIjkn  are 
constants.  Moreover,  we  take  t\  =  £2  =  1;  hence  Y  =  [0,1]  x  [0,1]  and 
|V'j  =  1,  |K'|  =  p (2  —  p)  so  that  0  ~  p(2  —  /i).  It  follows  that 

0 

- >2  as  p  — >  0. 

P 

One  can  establish  the  following. 


Theorem  4.2.  For  the  constant  coefficient  problem  with  =  i2  =  1,  there 
exists  a  sequence  pn  — ►  0  such  that: 

(i)  The  coefficients  p^jlh  °f  (\.o)  converge  in  the  sense  that 


/in  _  * 

P ij kh  *  Pi jkh 


where 

(4.8) 

Pijkh(s ) 


{ojjee  +  sbij({)(a((kh  +  sbakh) 
aaa  +  sl>wu 


(ii)  The  solutions  uMn  of  (HE  Mn  J  satisfy 


— >  u*  in  L°°(0,  T]  Hi (f!))  weak * 


where  u *  is  the  solution  to 

(HE) 

(2 phuwtt{t),ih)  +  <7*(t)(iT(-),</0  =  (2/,  0)  <  €  (0,T)  and  a//  0  €  //62(f>), 

u*(0)  =  lim  ^u0  =  u0/2 

u*(0)  =  lim  =  n0/2, 

M-*0  (/ 

with 

(4.9)  am{t)(um(-),if)=  f  [  Pijkhi1  -  T)u;kh{T,x)i!>,ij(x)dT  d.r. 

Ju  Jo  >oX,h 
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We  return  to  the  AFAL  grid  structure  of  Section  3  corresponding  to  a  di¬ 
rect  physical  model  using  the  Love-Kirchhoff  theory  with  Kelvin- Voigt  damp¬ 
ing  which  is  given  by  equation  (3.7)  with  the  sesquilinear  forms  given  by  (3.(i). 
This  was  the  abstract  variational  form  of  (3.1)  -  (3.5).  We  specialize  the  Je¬ 
suits  of  this  section  to  that  example.  We  have  the  coefficients  «7 "jt.'i  = 
given  by 

El 

ami  =  a2222  =  - - r  =  9l 

1  —  vi 

El 

au22  —  «22ii  —  v- - -  =  qiq  3 

1  —  vz 


a  1 2 1 2  —  a2i2i  —  ai22i  —  02112 


El 

2(l+e) 


with  all  other  atJkh  zero.  The  bl}kh  are  given  bv 

t  L  CD^ 

ami  —  02222  —  - - 7  —  <72 

1  —  vz 


5?i(l  -  +  ) 


hi  122  —  bj2i 1  —  ^ 


CIJ  I 
1  -  1+ 


=  9293 


bv. 


1212 


/2121 


—  hi  22 1  —  h/n-2  — 


CD  I 


2112  — 


2(l+i/) 


=  ij'/A  1  - 


with  all  other  b^h  zero.  In  this  case  (3.6)  and  (1.1),  (1.2)  agree.  Moreover, 
using  (3.8)  of  Theorem  1.2,  we  find 


Pim(5)  -  P2222 ( ■'' )  —  El  +  scD l  —  q\{  I  -  ql)  +  ,s q 2 ( 1  -  qr>.) 

Pnn(s)  =  F'2i2i('s)  =  Pnni*)  =  /+22i(-s)  =  TT~+STTT,  =  9i( 1  ~‘/a  )+■•“/*•<  1  ~  V  < 

with  all  other  p’jkh  zero.  It  follows  that  a*  Irom  (1.9)  can  be  written  as  the 
sum  of  two  sesquilinear  forms  rr*  and  rrm,  given  by 

o'* ( <P,  +)  =  Ji}{(l >0  —  9a )(o,i  1  C’,11  +O.22  t’,22  )  +  ‘b/,(  1  —  9:,) O.i 2  '-'02  }</.r. 

We  observe  that  a also  satisfy  the  conditions  (A).(E)  and  f(  )  of  Section 
2  so  that  the  convergence  results  summarized  in  that  section  are  readily  ap¬ 
plicable  to  the  homogenized  model  in  this  case.  If  we  write  the  homogenized 


model  system  in  strong  form  analogous  to  (3.1)  -  (3.5),  we  obtain  (we  include 
the  viscous  damping  term  so  as  to  compare  with  (3.1)) 


,dV  du*  d2  Ml  d2Ml2 

ph  dt 2  dt  +  dx\  +2dxxdx2 


d2  Ml 

+  dx\ 


t  €  (o,  T),  x  <=  n. 


(4.11) 


u*(0)  =  no/2,  u*(0)  =  t>o/2 


where  — >  u0,  — >  v0  weakly  in  L-i{Q)  as  t  — >  0  as  stated  in  Theorem 

4.1  and  the  moments  A/,1,  A//,  Ml2  are  given  by 


Eld2  a *  cDl  d3u- 

~2"dx[  +  IT  dx2dt 

El  d2u'  cDI  d3u 
2  dx\  ^  2  dx\ dt 

ei  d2u •  cD i  d3u- 

1  T  u  dx\dx2  1  4 -  v  dx\dx2dt 

We  note  that  M „12  is  the  same  as  A/12  of  (3.2)  while  the  form  of  equation 
(4.11)  is  the  same  as  that  of  (3.1),  with  only  the  moments  being  different. 
Of  course  u~  must  also  satisfy  the  clamped  boundary  condition  (3.3)  along 
the  edge  x2  =  0  of  fl  and  the  free  boundary  conditions  (3.4)  (with  A/J,  Ml2) 
along  Xi  =  0  and  x:  =  a  and  (3.5)  (with  Ml,  Ml2)  along  x2  —  6. 


Ml 

(4.12)  Ml 

Ml2 
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5.  Comparison  with  Experimental  Results  for  the  AFAL 
Grid 

The  homogenized  model  (1.11)  -  (4.12)  for  the  AFAL  Grid  is.  of  course, 
a  type  of  approximation  to  the  physical  model  (3.1)  -  (3.5).  Unlike  many  ap¬ 
proximate  models,  we  cannot  increase  the  order  of  approximation  by  chang¬ 
ing  some  parameters  or  mesh  sizes.  Even  though  (HE)  is  the  result  of  a 
limiting  procedure  c,p  — ►  0,  the  periodicity  of  bays  and  the  thickness  of 
members  in  a  given  grid  are  fixed  and  hence  c  and  p.  are  fixed.  We  can  ex¬ 
pect  the  homogenized  model  to  be  a  better  approximation  for  structures  with 
more  bays  and  thinner  members  than  for  structures  with  a  small  number  of 
bays  and  relatively  thick  members. 

Before  using  the  model  (4.11)  -  (4.12)  for  estimation  and  control  prob¬ 
lems  for  the  AFAL  Grid,  we  performed  some  initial  experiments  to  test  how 
well  the  homogenized  model  described  quantitative  properties  of  vibrational 
characteristics  of  the  grid.  One  such  investigation  involved  vibrational  ex¬ 
periments  (the  grid  was  displaced  initially  from  equilibrium  and  free  release 
vibrations  were  observed;  the  data  consisted  of  accelerometer  measurements 
at  several  locations  on  the  grid).  Analysis  of  the  data  included  experimental 
determination  of  the  fundamental  frequencies  (which,  of  course,  depend  on 
the  internal  damping)  for  the  grid.  We  then  compared  these  values  with 
those  predicted  by  the  analysis  of  the  homogenization  model  (4.1 1)  -  (4.12). 

The  AFAL  Grid  used  in  our  experiments  was  a  5 ft.  square  (a  =  b  = 
5)  with  16  =  4  x  4  square  bays  measuring  12"  x  12".  This  corresponds  to 
a  calculated  y  ~  1.2"  with  cf,  ~  14.4".  The  grid  was  constructed  from 
aluminum  alloy  6061-  16  for  which  hand  book  values  for  stiffness,  etc.,  were 
available.  We  used  ph  =  .0550132.s/u</.s///2,  v  =  5,  El  —  135.64 U>-  Jt.cp  — 
1.728  x  10°/6 •  se c/ft2  and  7  =  .02 slugs/ ft2  sec  in  our  calculations  with  the 
model  (4.11)  -  (4.12).  We  calculated  the  iirst  eight  frequencies  y  for  the 
homogenized  model  using  two  different  approaches:  an  approximate  mode 
shape  technique  based  on  textbook  approximations  [G],  [L]  and  an  eigenvalue 
analysis  (MATLAB)  for  a  finite  element  (bicubic  B-splines  modified  to  satisfy 
the  essential  boundary  conditions)  Galerkin  approximation  to  (4.1  1)  -  (4.12). 


1!) 


We  obtained  essentially  the  same  results  with  both  approaches.  These  values 
are  listed  along  with  the  experimentally  obtained  frequencies  u>fxl  in  the 
table  below.  We  observe  that  for  i  =  1,2,3,  there  is  good  agreement  while 
ujI  is  approximately  the  average  of  u'J  xp  and  ivPXP.  There  is  reasonably  good 
agreement  between  the  pairs  u*  and  u^xp  for  i  =  5,6,7.  This  comparison 
suggests  that  use  of  the  homogenized  model  is  a  reasonable  approximation 
at  least  for  vibrations  involving  modes  with  frequencies  less  than  16/ic. 


le  i 

< 

,EXP 

u,, 

1 

.785 

.781 

2 

2.52 

2.15 

3 

4.93 

4.69 

1 

6.84 

6.35 

5 

8.44 

7.23 

6 

13.8 

8.11 

7 

15.1 

13.67 

8 

15.5 

15.5 

Table  4.1 
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